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Abstract: 

In this article, the existence of gravity-capillary waves travelling down the surface of a falling 
liquid film past a vertical plate has been considered. Kapitza's scheme of finding an 
approximate expression of velocity u when the film surface assumes an arbitrary shape y = h(x, 
t), which changes with time, has been emphasized. The expressions for dimensionless 
wavelength, dimensionless wave number and Weber number have been obtained and are 
computed for an admissible range of the wave celerity. The stream line pattern has also been 
studied and presented through graphs. 

Key words: Wave motion, liquid film, gravity-capillary waves, stream line pattern 



Department of Mathematics, University of Botswana, Private Bag 0022, Gaborone, 
BOTSWANA 

A Quarterly Double-Blind Peer Reviewed Refereed Ope n Access Internati onal e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., W HJiMdMJIIiMM as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Engineering, Science and Mathematics 
http://www.ijmra.us 





Volume 2, Issue 3 



ISSN: 2320-0294 



INTRODUCTOION 



The flow of a thin liquid is a special case of flow in which both surface tension and 
viscosity play dominant role and is of fundamental importance in distillation of gases and 
liquids, having varied applications in engineering operations such as distillation, gas absorption, 
condensation of vapors, where heat and mass transfer are intimately connected to fluid motion. 
Extensive reviews by Norman [1], Dukler and Wicks [2], Levich [3] and Fulford [4] bring out 
the importance of such studies. 



The experimental studies of the flow of thin films have revealed that (i) at Reynolds 

numbers Re = U °^° that do not exceed 20 to 30 (where u is the average velocity over the cross- 
v 

section of the film of thickness h and v the kinematic coefficient of viscosity), there exists the 
usual viscous flow regime and the film thickness is constant, (ii) at Reynolds numbers Re > 30 to 
50, a so-called wave regime appears, in which the wave motion is superposed in the forward 
motion of the film and (iii) at Reynolds numbers Re about 1500, the laminar regime breaks down 
and the turbulent motion sets in. In the case of very thin films, the break-up of the film into 
separate drops is caused by the action of capillary forces. The direction depends on wetting 
conditions. 



The undulatory character of the flow was experimentally detected by Kirkbridge [5]. 
Their studies reveal that this motion appears for values of Reynolds number about 20 to 30. 
These studies show that the law for laminar flow is obeyed only as an average and refers to the 
mean thickness of the layer. In fact, the character of the flow differs from that of a simple 
laminar flow. The reason for this difference is sought due to the fact that in deducing the Navier- 
Stokes equations, one does not take into account surface tension which in a flow of a fluid with a 
low viscosity and in thin layers, even if the free surface is slightly deformed, acquires a 
considerable value quite comparable with forces of viscosity. If capillary forces are taken into 
account, the undulatory flow even at small velocities is indeed more stable than a simple laminar 
flow. The existence of a more stable undulatory flow is of great interest because it allows one to 
explain and describe a number of known physical phenomena observed in the flow of thin layers. 

The theoretical aspects of these flows fall mainly under three categories, namely, (i) the 
search for periodic steady states (ii) the stability analysis and (iii) the semi-empirical 
characterizations. Kapitza [6] was the first to predict the appearance of gravity-capillary wave 
theoretically and the results of theoretical analysis of the wave regime in the flow of film down a 
plate were verified, experimentally, by Kapitza and Kapitza [7]. Direct measurements of wave 
profile made by the shadow graph optical technique confirmed the theoretical predictions, 
particularly that of the amplitude. It was found that wave motion of the film readily gives way to 
turbulent motion when there are external disturbances to the liquid flow. 

Chaturvedi et al. [8] have considered exponential fluid flow along an infinite porous plate 
with constant suction. They observed that there is no back flow near the wall. Raju [9] has 
studied the oscillatory flow past an infinite vertical porous plate. He revealed that in the case of 
heating of the plate, the mean velocity increases whereas it decreases in the case of cooling of the 
plate with increase in Prandtl number. Masthanrao et al. [10] have studied chemical reaction and 
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hall effects on MHD convective flow along an infinite vertical porous plate with variable suction 
and heat absorption. They observed that an increase in the Prandtl number leads to a decrease in 
the primary and secondary velocities, and also a decrease in the primary and secondary 
temperatures. Manjulatha et al. [11] have studied the effect of aligned magnetic filed on free 
convective flow past an infinite vertical flat plate through porous medium with temperature 
dependent heat source. They found that the temperature and velocity increase with the increase 
of heat source parameter. 

One of the earliest contributions to the analysis of a falling liquid film (Fig. 1) was made 
by Nusselt [12]. The equation for slow steady state of a film under the action of gravity is 
given by 



d u _ 
dy 



(1) 



The boundary conditions are: 



(i) at the free surface of the liquid, y = h , ju — = , and 

dy 

(ii) at the wall y = 0, u = 0. 



y = h (x,t) 




Fig.l Physical model 
The solution of (1) under the above boundary conditions is 
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3u 



y 



K V 2K 



(2) 



where u a = — \udy = ^— £ - (3) 

h o * 3v 

This theory is invalid when waves appear at the free surface. But it is extensively used as a first 
approximation. 

Kapitza [13] made an attempt to describe the wavy flow by postulating an oscillatory steady state 
solution of the equations of motion. His theory was described in detail by Levich [3] and 
Bushmanov [14]. All these researchers considered the hydrodynamic equations of motion and 
boundary conditions for thin liquid film and took into account the effects of surface tension. 

However, they have neglected the term v as in boundary-layer theory and variation of the 

dx 

surface coordinate h with the longitudinal coordinate x at several stages in their derivation of the 
differential equation, Levich [3] : 

/ch a d 3 y/ , .. 9 .dy/ 3v , . . 3vu a n 

— —rY + (c-u )(c-—u )^-—(c-3u )y/+g — -^- = (4) 
p dx 10 dx h o h o 

where 

h = h a (l + y/) (5) 
and \\f the steady state periodic solution of the equation (4) is given by 
In 

y/ = Asin — (x-ct). (6) 

They argued that in order that the equation (4) to have undamped periodic solution, both the free 
term and the coefficient of \\i must be equal to zero. 

The first of these requirements is that 

«-^ = (7) 
h 

which signifies that in a first approximation, the thickness h a is close to the thickness of the layer 
in laminar flow down the plate. 

The second requirement takes the form 
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c-3w o =0 (8) 
and defines the phase velocity of the wave in equation (6) as the first approximation. 
The frequency of the wave w is given by 



v = — c = 3u 



4.2 p 

Kh„ 



(9) 



An important feature of this equation is W = U and that both approach zero simultaneously. 
This is the basic difference between surface waves on a film flowing down a plate and ordinary 
capillary waves. The undamped nature of waves on a film of viscous liquid flowing down the 
plate also differentiates them from the ususal capillary waves. Waves on the film surface are 
constantly maintained by the force of gravity, which moves the film as a whole. Kapitza, next 
obtained a second order of approximation, a non-linear treatment but lacks justification 
according to Massot et al. [15]. A quite different periodic solution was obtained by Ishihara et 
al. [16] who solved the hydraulic equation of motion for unsteady flow in open channel with a 
constant inclination. Neglecting the surface tension effects, they obtained 



c — u 



6 - + 

5 V 25 mRe 



(10) 



where m is the slope of the channel. 



Massot et al. [15] elaborated the work of Kapitza by considering the complete x-component 

equation of motion without neglecting ^-^ , the complete y-component equation of motion and 

dx 

took into account the variation of h with x, explicitly at all stages. Their results thus include as a 
special case of Kapitza and Kapitza [7] and of Ishihara et al. [16] in the frame work of linear 
approximation. 

Yih [17], Benjamin [18] and, Jones and Whitaker [19] investigated the stability of 
Nusselt's solution in terms of the linearized theory of small perturbations. By postulating that 
waves are generated by infinitesimal periodic disturbances, and that the most likely to be seen 
are the fastest growing ones, this theory can be utilised to obtain a partial, if not satisfactory, 
description of the wave motion. 

Philips [20] has examined wavy films from a statistical stand point and Lee [21] used the 
semi-empirical turbulent model developed for flow in closed conduits (Lee, [21]). However, 
these works are not detailed further here, since we are considering the apparently laminar range 
in which regular wave occurs. 
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2. THE MATHEMATICAL METHOD 

2.1 The basic equations 

The development of this method starts with the equations of continuity and motion for 
two-dimensional vertical flow of an incompressible Newtonian liquid of constant viscosity and 
is based on the search for an oscillatory steady state solution. 



The equation of continuity is 

du du 
— + — = 0. 

dx dy 

The x-component of Navier-Stokes equation is 

du du du ( d 2 u d 2 u 



1 dp 
+ g — . 

dx 2 dy 2 ) p dx 



(11) 



(12) 



dt dx dy 

We propose to solve these equations (11) and (12) subject to the boundary conditions 

— = at y = h(x, t) for all x and t, (13) 

dy 

u = v = 0aty = for all x and t, (14) 

and 



d 2 h . x 
P = Pen ~ K — j at v = n(x, t) 
dx~ 



(15) 



where k is the coefficient of surface tension. The radius of the curvature of the film is 

d 2 h 

approximated by , since the surface slope is expected to be so small that its cube can be 
neglected. 

In addition to the equations (11) and (12), we add the boundary-layer approximation of 
hydrostatic equilibrium across the film to be valid, following Kapitza [13] and Levich [3]: 



^ = 0. 

dy 



i.e., p = p(x, h(x, t), t). 



(16) 



As a matter of fact, this is a consequence of the y-component of Navier-Stokes equations of 
motion Schlichting, [22]. 
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dp 

dx 



572 



With this, the equation (12) becomes 

du du du ( d 2 u d 2 u 



, rg+- 

dt dx dy y dx dy J p dx 
The surface coordinate h is related to u by the mass balance equation 



k d 3 h 



dh_ 
dt ' 



judy = 



-(uh) 



Sx- dx 

where u is the average velocity over a cross-section defined by 



u= ] -'\udy. 



(17) 



(18) 



(19) 



hi 

The equations (1 1), (17) and (18) are the general equations of motion of the liquid in a thin film. 
2.2 The derivation of the wavy equation and its compatibility 

We follow the Kapitza's scheme of finding an approximate expression of u when the film 
surface assumes an arbitrary shape y = h(x, t), which changes with time. We note that u must 
satisfy the conditions (13) and (14). These conditions led him to pattern the approximate 
expression after the Nusselt flow. In the same spirit, we write the following expression for u for 
the flow under consideration. 



3u(x,t) 
' Kx,t) 



2h(x,t) 



(20) 



We substitute (20) in the equation of continuity (11) and integrate under the boundary conditions 
(14) and v(x,0,t) = 0. Then 



(21) 



3 du\ 




3u(x,t) dh(x,t)\ 


V y 3 


h(x,t) dx\ 


_6h(x,t) 2 _ 


\ + h 2 (x,t) dx 


2 3h(x,t) 



Substituting (20) and (21) in the equation (17) and integrating the result with respect to y from 
to h(x, t), we get 
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du 1 u dh 9 - du 
dt 2hdt lo" dt 



3 u <3/z 
10 h dx ' 



^ k d 3 h + | 5 2 m 
/? dx 3 \ dx 2 



1 m d 2 h _ 3u _ 5 1 du dh) 
2hdx 2 h 2 Ahdxdx) 
(22) 

The principle here is to determine u so that equation (20) becomes the best representation of u in 
the sense that it satisfies (22) in an averaged manner. It may be recalled that in this respect, it is 
like the Kerman-Polhausen technique in solving boundary-layer flow problems, Schlichting [22]. 
If a solution of equation (22) exists such that it represents a steady travelling wave in the x- 
direction (i.e., periodic in x and t), it must be a function of a simple variable x - ct, where c is the 
phase velocity of the progressive undamped wave train. To search for such a solution, we can 
replace dl dt in equation (22) by -c(d/dx) and linearize the equation (22) on the assumption 
that the wavy surface is but a small perturbation on the flat surface y = h , i.e., 



h = h a [l + w(x,t)] 



(23) 



where < 1 . 

The mass balance equation (18) gives us 
c{h-h n ) = uh-uh 

where u is the average velocity of the liquid film in equation (3) corresponding to Nusselt 
solution. 

In terms of y/ , this gives us 

u = u a +(c-u )y/ . 

In a similar manner, 

du . dw 

— = -cic-uj— 1 —, 
dt " dx 



du . ^dy/ 
— = (c-u )— ! -, 
dx dx 



-[u o +(c-2u a )i//], 



h h 
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dh , dw 
— = -ch D — L - 
dt dx 



Substituting these equations in the equation (22), we arrive at a linear equation, which is 
expressed in a non-dimensional form using 

u h 

Re = - £ - £ - (Reynolds number), (24) 



Fr = — — (Froude number), (25) 

sK 



A, log 

- — J (dimensionless wavelength), (26) 

2n V k 



N = (dimensionless wave number) (27) 

We = = P 2 N 2 Fr (Weber number) (28) 

Kip 

(29) 



and 



A 

?] = — (Wave celerity) (30) 



where X is wavelength. Thus, we get 



^^ + ^7-|^V" + ^ 2 -y7 + |]^^^-(7-3V = (31) 

where the primes represent differentiation with respect to £ Equation (31), being homogeneous, 
will admit a sinusoidal solution, only if 

P'FrJrj'-^rj + fi 1 (32) 
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3(^7-3) 
3 ' 

v-- 



P in equation (32) is real, if 

rj > T] 2 or rj<rji 
where r|i < r\2 are the roots of the quadratic equation 



2 12 6 _ 

77 ri + - = 

5 5 



(33) 



(34) 



ft- 1— J and ^H— ■ 



On the other hand iV in equation (33) is real, only if 
3 

— <T]<3. 



(35) 



Thus an undamped, steadily travelling wave in the form of a sinusoidal function is compatible 
with the physical situation if positive values of r\ lie within the common range as indicated in 
equations (34) and (35). No solution in the particular form under consideration exists if the two 
regions do not overlap. 



It is clear from equations (34) and (35), that the common range for r\ is 
1.689 <?7<3. 

For each such admissible r\, we can find the wavelength of the corresponding wavy flows 
277-3 * 



6 + V6 
5 



X = 2nh ( 
and the Weber number is 

We = - 



6(77-3) 



30(3-77) 



~(2t7-3)(5t7 2 -12?7 + 6)' 
3. THE DESCRIPTION OF THE WAVY FLOW 



(36) 



(37) 



A Quarterly Double-Blind Peer Reviewed Refereed Ope n Access Internati onal e-Journal - Included in the International Serial Directories 
Indexed & Listed at: Ulrich's Periodicals Directory ©, U.S.A., W HJiMcMdMlliMM as well as in Cabell's Directories of Publishing Opportunities, U.S.A. 

International Journal of Engineering, Science and Mathematics 
http://www.ijmra.us 




Volume 2, Issue 3 



ISSN: 2320-0294 



3.1 Dimensionless wave celerity 



At low Weber number (thin films), we have 77 = 3 which is the value predicted by the 
theories of Kapitza [13], Yih [17], Benjamin [18], Hanratty and Hershman [23]. Thus, these 
theories are asymptotic to the present solution at zero Weber number. 



Also when 77 — > 1.689, We — > 00. This result is the limiting value of Ishihara et al. [16] 
for large Weber number (thick films). Their analysis is therefore asymptotic to the present 
solution at high Weber number. 

The solution in equation (37) gives the variation of the Weber number We with the 
dimensionless celerity?;. Fig. 2 represents a comparison of all the existing theories. Further, at 
medium Weber number, where the other theories fail, there is a qualitative agreement between 
experimental results and the present predictions. 




Fig. 2 Weber number We as a function of wave celerity 77 
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It is remarkable that all experimental observations of 77 lie between 3 and 1.689, which is the 
predicted range. Table 1 gives the computed values of We and N for the admissible range of r/. 
The values of Weber number We and wave number N vary from zero to infinity as the wave 
celerity 77 varies from 3 to 1.689. 

TABLE 1 

Computed values of Weber number dimensionless wave number for admissible values of 
wave celerity 



3.0 
2.9 
2.8 
2.7 
2.6 
2.5 
2.4 
2.3 
2.2 
2.1 
2.0 
1.9 
1.8 
1.7 



0.0000 
0.0809 
0.1989 
0.3731 
0.6342 
1.0345 
1.6666 
2.7062 
4.5113 
7.8948 
15.0001 
33.0006 
100.0059 
1950.7978 






0000 





4269 





6793 





8660 


1 


0445 


1 


2248 


1 


4142 


1 


6202 


1 


8516 


2 


1213 


2 


4495 


2 


8723 


3 


4641 


4 


4159 




00 



3.2 Dimension wavelength 

The dimensionless wavelength P comes out as a function of the Weber number through 
the Reynolds number 



Re = ^'{^T 



(38) 



for assumed values of fluid parameters such as p, k and //. 
From equation (32), we have 
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when r] — > 3, We — > 0, P — > oo and when 77 — > 772 = 1 -689, We — > 00, P — > 00. 

The equation (39) is used to plot the dimensionless wavelength P against the Weber number We 
in Fig. 3 for alcohol and water. 




id 3 10" 2 10" 1 1 10 id 2 10 3 10 4 10 s 10 6 



Fig. 3 Dimensionless wave length P as a function of Weber number We 

It is observed that the wavelength passes by a minimum (this minimum occurs when 77 = 1.9 for 
any liquid this value corresponds to the Weber number 33 and the dimensionless wavelength P = 
0.38 (for water) and P = 0.59 (for alcohol) and becomes very large when We — > and We — » 00. 
Such a behaviour has been observed experimentally by Kapitza and Kapitza [7] whose 
experimental results seem consistently higher than the theoretical predictions. 

3.3 Dimensionless wave number 

In a similar way, the dimensionless wave number /V is plotted against the dimensionless 
wavelength P in Fig. 4. The path of the system for increasing Weber number goes through three 
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regions related to Robertson [24] to different known types of waves namely: long waves, 
gravity-capillary waves and deepwater waves. The transition between the gravity and capillary 
waves is usually set in the literature at P = 1, because this value corresponds to the minimum of 
celerity of an ideal liquid. 




Fig. 4 Dimensionless wave number N as a function of dimensionless wavelength p 



4. STREAM LINE PATTERN 

The stream lines are defined by 

dx dy 



(40) 



Let S = y/h denote the dimensionless transverse coordinate for stream line calculation. At a 
fixed instant t (t = say), using the equations (20) and (21) for u and v in (40), we obtain 



25' (3-S)' 
8 + 3-S 



(uh)' 
nil 



(41) 
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where the primes represent differentiation with respect to x. On integration, we find that the 
dimensionless stream function *F (with respect to u h ) is 



*P = (3- S)\l + ?7Asin^ 



(42) 



in which the form of y/'m equation (23) is 
y/ = Asm . 



(43) 



Two families of stream lines appear depending upon the dimensionless celerity - 
amplitude number C = rjA being greater or less than unity. When C < 1, we obtain the 
unidirectional stream lines. Fig. 5 gives a general view of the instantaneous stream lines of a 
falling liquid film corresponding to a slightly deformed parallel flow and to positive values of the 
stream function y/. 



Fig 5. Streamline pattern of a falling liquid film for C < 1 showing a unidirectional flow 
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When C > 1, we obtain reversal flow stream lines, a general view of which is shown in Fig. 7. 

Fig. 8 gives a detailed stream line when y/= -2.4, -1.0, 0.0, 1.0, 5.0 and 8.4 for C = 1.8. Thus 
when C > 1, we obtain reversal stream lines corresponding to positive and negative values of the 
stream function, suggesting a swirl in the liquid. 




Fig. 8 Streamlines for C= 1.8 >1 

Using equations (18), (19) and (23), we obtain 

- cu/ + u 1 

u =— -. 

l + y/ 



(44) 



The transition between the above two types of flow is determined by setting u = . We obtain, 
on using equation (43) 



Csin^ = -1. 
1 



(45) 



These two types of flow may be expected to reveal different mechanisms of surface renewal. 
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5. CONCLUSIONS 



We have observed that for thin films (i.e., at low Weber number), we have the wave 
celerity rj = 3 which is the value predicted by the theories of Kapitza, Yih, Benjamin, Hanratty 
and Hershman and hence these theories are asymptotic to the present solution at zero Weber 
number. Also, it has been found that when 77 — > 1.689, We — > 00 which is the limiting value of 
Ishihara et al. [16] for thick films (large Weber number). Their analysis is therefore asymptotic 
to the present solution at high Weber number. At medium Weber number, where the other 
theories fail, there is a qualitative agreement between experimental results and the present 
predictions. 



It is remarkable to note, from Fig. 2, that all experimental observations of rj lie between 3 and 
1 .689 and this is the predicted range. It is noticed from Table 1 that the values of Weber number 
and wave number vary from zero to infinity as the wave celerity varies from 3 to 1.689. From 
Fig. 3, it is observed that the wavelength passes by a minimum when 77 = 1.9 and becomes very 
large when We — > and We — > 00. Such a behavior has been observed experimentally by Kapitza 
whose experimental results seem consistently higher than the theoretical predictions. It has been 
found from Fig. 4 that the path of the system for increasing Weber number goes through three 
regions to different known types of waves namely: long waves, gravity-capillary waves and 
deepwater waves. It is found that from Figs. 5 and 6, when C < 1, we obtain the unidirectional 
stream lines. When C > 1, we obtain reversal flow stream lines, a general view of which is 
shown in Fig. 7. Also from Fig. 8, it has been observed that when C > 1, we obtain reversal 
stream lines corresponding to positive and negative values of the stream function, suggesting a 
swirl in the liquid. 
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